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Abstract 

Chiral effective field theory predicts a specific charge symmetry violating ampli- 
tude for pion production. This term is shown to provide the dominant contribution 
to the forward-backward asymmetry in the angular distribution for the reaction 
pn — > dtr , for reasonable values of the mass difference between down and up quarks, 
Smjy- Using a value 5m n ~ 3 MeV leads to a prediction of a 10 standard deviation 
effect for a TRIUMF experiment. 
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Charge symmetry is an approximate symmetry of QCD. If one ignores electromagnetic 
corrections and the mass difference between u and d quarks, then the Lagrangian is 
invariant under the interchange of u <-> d. Isospin invariance implies charge symmetry, 
but the converse is not true. 

Much has been learned about charge symmetry breaking (CSB); see the reviews |], |f. 
The general importance of the quark mass difference effect was stressed in Ref. J2|. Our 
purpose here is to use chiral effective field theory to make predictions, based on QCD, for 
the reaction np — > dir°. 

The use of chiral Lagrangians, which employ hadronic degrees of freedom, represents 
a serious effort to provide a rigorous and model-independent methodology to use QCD 
to make predictions at relatively low energies 0, ||, This is because the most general 
Lagrangian which respects unitarity, has correct properties under cluster decomposition 
and the same symmetries as QCD, should give results that encompass the results from 
QCD proper ||. Predictive power is retained because at low energies an expansion in 
momentum can be formulated using power counting arguments. 

We wish to exploit the feature that chiral symmetry makes predictions for effects that 
stem from the small but explicit breaking of chiral symmetry generated by the quark 
masses. Consider the proton-neutron mass difference m n — m p , which includes a contribu- 
tion Smjsr from the quark mass difference — m u = e(mrf + m u ) (e ~ 1/3), and another 
contribution Sm^ from electromagnetic effects. Thus m n — m p = 5m^ + 5m^. Neglecting 
strangeness, chiral symmetry is essentially an SO (4) internal symmetry. One can show 
that the quark-mass-difference term in the QCD Lagrangian behaves under SO (4) as the 
third component of an 50(4) vector ]7], || ||. Therefore, in the effective hadronic theory, 
all isospin- violating interactions generated by the quark mass difference break 50(4) as 
third components of (tensor products of) 50(4) vectors. The operators of interest here 
involve the nucleon N and pion ir fields. The leading isospin-violating term at low en- 
ergies coming from the quark mass difference is the third component of an 50(4) vector 

11 = 

£fi = ^ (nW 3 N - ^JVWr ■ riv) , (1) 

where r represents the Pauli matrices in isospin space, F n = 186 MeV is the pion decay 
constant, and D = l+iz 2 /F%. We use D — 1 here. Eq. (P represents the isospin-violating 
part of the nucleon sigma term. 

Furthermore, one can show that quark interactions generated by ("hard") photon ex- 
change break 50(4) as the 34 component of an 50(4) antisymmetric rank-2 tensor [Q. 
Thus the low-energy effective theory must include isospin-violating interactions that break 
50(4) as 34 components of (tensor products of) 50(4) antisymmetric rank-2 tensors. The 
leading isospin-violating term of this kind at low energies was shown to be ||: 

4 P 1} = ^ (Vr 3 iV + • r - 7rV 3 )iv) . (2) 

Observe that chiral symmetry links the first terms in Eqs. (0,§), which are contri- 
butions to the nucleon mass difference, to the second terms, which are isospin-violating 
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pion-nucleon interactions. Since the parameter F w is determined from another process 
(pion decay), these pion-nucleon interactions must exist with the given strength if our 
understanding of QCD is not totally flawed. 

Let us discuss the values of the parameters, 5m N and Sm^, which are not determined 
by chiral symmetry. The sum of these terms is fixed by the nucleon mass difference 
m n — m p = 1.3 MeV. Dimensional analysis suggests their separate order of magnitude in 
terms of the natural QCD mass scale, Mqcd ~ 1 GeV. We know 5rriN oc rrid — m u and 
oc m d + m u , so 5m N = O^em^/ Mq C d) ~ 7 MeV. The electromagnetic term Stun, 
long known to be negative ||10|| , should be proportional to the fine structure constant a, 
so that Sm N = —O^aMqcD/^) ~ —2 MeV. Estimates such as these cannot be trusted to 
better than a factor of a few. More precise estimates can only be made at present with 
model-dependent assumptions. Most models estimate that Sttln ~ —a/R,N — —1.5 MeV 
0, in which case 5rriN — 3 MeV. We will refer to these values, which are also consistent 
with the results of Gasser and Leutwyler [JTTJ] , as "quark model estimates" . 

Verifying the existence of the terms of Eqs. (|l],|2|) is of high interest, as is determining 
the parameters. The simplest procedure would be to study pion-nucleon scattering close 
to threshold [|7|, |8], |9|, [12| . The lack of 7r° beams requires the comparison of charge exchange 
and elastic scattering processes. If one neglects electromagnetic effects and the kinematic 
dependence on the threshold values of the energy on the masses, the leading-order chiral 
Lagrangian gives for the triangle discrepancy among the pion-nucleon amplitudes T: 

^T(tt + p — > 7i + p) — T(n~p — > ix~p) — \/2T(ix~p — > n°n) ^ Sm N — Sm N 
T(n+p -> n+p) - T(ir-p -> z~p) + V2T(ix-p -> ix°n) ~ 2m w 



The quantity D is found [[13] to be about 6%, leading to a very large negative value of 



drriN — drriN ~ — 17 MeV. If the 6% result is correct, there seems to be a serious problem 
with our current understanding of low-energy QCD. However, electromagnetic effects 
are very important at low energies as are the kinematic difference in threshold energies. 
Kinematic effects and sub-leading strong interactions have been studied in Ref. [fL2|, but 



a precise accounting for the electromagnetic interactions presents a very serious challenge 
to theorists. A small shift in the normalization could be reflected in a huge change in D. 
It seems that a successful determination of Stun and Em^ from ttN scattering will require 
a strong theoretical and experimental effort. 

The biggest problem in all of this is the lack of a pion beam. The reaction np — > dir 
is relevant because it involves the production of a virtual pion from one nucleon ("pion 
bremsstrahlung" ) , followed by rescattering on the second. Charge symmetry predicts that 
the angular distribution is symmetric about 90° in the center-of-mass. Any asymmetry 
must be caused by charge-symmetry-violating effects. We here examine the integrated 
forward-backward asymmetry in the center of mass of the reaction np — > dn , 

ff da \a(0) - <r(x - ff)\ 

JpfM«) ■ (4) 

Our interest in this particular observable is motivated by an experiment in progress at 
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TRIUMF [14] which aims to measure this asymmetry close to threshold, = 279.5 
MeV or 77 = q^/m^ = 0.17, with a precision of ±0.06% []. 

A arises from an interference between s- and p-wave pions and thus decreases as 
one approaches threshold energies. Nevertheless, there are advantages to considering 
this pion production reaction at threshold kinematics. Near threshold one minimizes the 
four-momentum transferred k, enhancing isospin violation relative to chiral invariance, 
because chiral (thus isospin) conserving contributions to ttN rescattering vanish with 
the pions' four-momenta. Right at threshold the leading ttN rescattering mechanism 
( Weinberg- Tomozawa term goes as the sum of the pions' energies, uo q + u^, the real 
pion has as small energy as possible, u q ~ m^, while the typical energy of the virtual 
pion is just Uk ~ m n /2. In pion production the three-momentum p w y/m^m^ of the 
virtual pion is relatively large, which slows the convergence of a momentum expansion. 
However, there is a trade-off here, since the large momentum reduces the importance of 
electromagnetic effects. Note, also, that at any given order in a momentum expansion 
in p/Mqcd, the number of relevant isospin-breaking contributions is more limited than 
isospin-conserving ones. This leads to an isolation of Stun and 5m n effects. Even if other 
mechanisms (e.g. tc-t] mixing) are also relevant, they are not expected to exactly cancel 
among themselves. Thus, a model adjusted to correctly describe the total cross-section 
might provide a reliable estimate of the overall size of the 5m jy and 5m n contributions to 
the asymmetry; and these estimates would be confronted with data. 

The most recent calculation of the near-threshold asymmetry |jl6 includes effects of 



the nucleon-mass splitting on the ttN vertex function, which influences both the pion 
production amplitude and the nucleon-nucleon potential, and the far larger influence of n- 
•q-rf mixing, both in the pion-production kernel and in initial and final interactions between 
nucleons. This calculation does not include the seagull interactions ([JJ^) proportional 
to 5m N and 5m N . The isospin-conserving interactions are included in a nucleon-delta 
coupled-channel model that describes both pion production and iViV phase shifts to an 



accuracy of a few degrees from threshold over the delta region [17]. Close to threshold 
77-77 mixing is the dominant effect, taking place both in pion s and p waves. 

Here we estimate, using both a power counting argument and a numerical evaluation, 
the effect of 5m^ and 5m^ in the asymmetry (^) close to threshold. These mechanisms 
appear in the s wave, and can be relatively large because the isospin-conserving p wave 
amplitude is not vanishingly small when compared to the isospin-conserving s wave, even 
below 300 MeV. A power counting argument suggests immediately that the contributions 
from Eqs. ([DID should be dominant. Let us take the contribution from the 5m^ seagull, 
Fig. [I|a, as an example. The isospin-conserving p wave emission is proportional to q/F n , 
while the isospin-violating s wave gets a p/F n from the virtual pion emission, 1/p 2 from 
the pion propagator, and Sm^/F 2 from the seagull vertex. The product of the two 
amplitudes is then expected to be 

m n \ 3 / 2 q 



l The definition of the asymmetry Afi, used here is 1/2 of the definition used in Ref. p3 
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Figure 1: Contributions to A^: (a) -M qm ; (b) M.^. Here nucleon (pion) propagation 
is represented by a solid (dashed) line; and pion emission/absorption is denoted by a 
(circled) solid dot in case of an isospin-conserving p (s) wave, and a (circled) cross in case 
of an isospin- violating s (p) wave. 



in which the notation qm is used to denote the effects of Eq. ([[]). On the other hand, the 
contribution from 7t-i] mixing is smaller by a factor of m n /mN near threshold. This can be 
seen by integrating out the rj. This produces an isospin-violating pion-nucleon coupling 
of size /3q/F n in the p wave and (3ui q p/ F w m N in the s wave, where j3 ~ em^./m 2 N |19| 



Let us consider now isospin violation in the p-wave pion emission, Fig. [T]b. The isospin- 
conserving s-wave emission is proportional to uo q p/ 'F^m^ , while the isospin-violating p 
wave gets a p/F n from the virtual pion emission, 1/p 2 from the pion propagator, another 
p/Fn from the virtual pion absorption, l/u q from the nucleon propagator, and f^q/F^ 
from the real pion emission. The product of the two amplitudes is then expected to be 



\ 5 / 2 q 



The same estimate holds for the interference between isospin-conserving p wave and 
isospin-violating s wave. A comparison of Eqs. @ and @ reveals the ratio m n /mN 
between the two effects. Note that the counting rules used here are consistent with those 
of previous works in pion production pi 



We calculate A^ following the procedure described in Ref. |L6[ and earlier references. 
Before presenting the results, we display aspects of the relevant operators. Denoting the 
two nucleons with superscripts W and ^ , the spin-isospin structure of the isospin-violating 
rescattering is (at threshold, and neglecting terms of order m n /miy ) 

/ 



O qm = - - I ^ [SmNir^ ■ t< 2 > + r 3 (1) r 3 (2) ) - 8m N (r^ ■ - r 3 (1) rf)] 



d 



-fir 



where \x is the charged pion mass, p! = V3p/2, and / is the pion-nucleon coupling constant 
(/ = gAp/F n in leading order; / 2 /47r = 0.075 is used here). The isospin operator allows 
transitions from states with isospin / = to the deuteron, also 1 = 0: 

(I = o| T « ■ ± tPtP\I = 0) = -3 T 1, (8) 
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where only the isospin is indicated. Thus there is a large matrix element. If we do the 
spin-angle evaluation of the matrix element of O qm between 1 P 1 and 3 Si (for the deuteron) 
states we find 



-fi r 



We emphasize that the asymmetry in pion production thus depends on a different 
combination of 5mN and Sm^ than either the total nucleon mass splitting or the isospin 
violation related to charge-exchange pion-nucleon scattering. Under the assumption that 
we have the isospin-conserving mechanisms of pion production under control, can be 
used in conjunction with m n — m p to extract Sm^ and 5m^. 

The above operator may be compared with the operator involving phenomenologi- 
cal isospin-conserving, on-shell pion rescattering, which in leading order is given by the 
Weinberg- Tomozawa term []15|] , sometimes written in terms of A2 = /i 2 /47rF^ « 0.045. 
If we consider threshold kinematics, and only the matrix element between 3 Pi and 3 Si 
states, we obtain |20| 



°" = ~2^m Xi !(T< " x T<2))3(<r<1) + " m) ■ 1 — 1 • (10) 




(■wiPfli^^^p). (u: 



so that 



Assuming that the initial np 3 Pi and x Pi wave functions are the same (the phase shifts 
differ by only ~ 10% at our energies) so that the radial matrix elements of the same, we 
obtain a rough estimate for the ratio of matrix elements: 



fSilO^Pi) 2V2Sm N 



( 3 Si\O ic \ 3 Pi) 3 



0.03, (12) 



using the quark model estimate Sm^ — 6171^/2 ~ 4 MeV. This comes from taking Sm^ = 
3.1 MeV@ so that Sm^ = —1.8 MeV. This provides a first estimate of in the percent 
range. (Cf. Eq. (H) for the triangle discrepancy.) 

We turn now to the results of our detailed numerical evaluations. For illustration here 
we again use the quark model estimate 5m jv — 5m^/2 ~ 4 MeV. The amplitudes arising 
from using this value in Eq. (^) (denoted as loosely as "qm"), as well as those obtained 
from the sum of the previously computed amplitudes of Ref . |L6| (denoted loosely as ^77) , 
are listed in Table These are reduced matrix elements in the definition used e.g. in 



Ref. ||21|| . Some overall kinematic factors are ignored. 

At 279.5 MeV, the energy of the TRIUMF experiment, the computed contribution from 
5m N and Smjy to the asymmetry of Eq. (|) is found to be 0.97%, or about -3.5 times the 



strength of (and opposite in sign to) the sum of the contributions calculated in Ref. |16 | 



This earlier calculation gave Af b =—0.28%. As a consequence, the total asymmetry is 
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Table 1: Real and imaginary parts of amplitudes at rj = 0.17, E = 279.5 MeV. The CSB 
amplitudes have been multiplied by a factor of 10 3 . See the text for the definition of qm 
and r/iT. 





qm amplitude 


7T7] amplitude 


transition 


real part 


imag part 


real part 


imag part 


1 P l -> s 


25.06 


-13.25 


-10.38 


6.34 


3 S 1 ^p 


3.01 


0.68 


-3.86 


-1.36 


3 A ^ P 


0.82 


-0.10 


-7.26 


2.65 


3 D 2 ^p 


6.47 


2.76 


2.89 


-0.3 




0.901 


-0.487 


0.901 


-0.487 


1 D 2 ^p 


0.535 


0.086 


0.535 


0.086 



Aji, = 0.69%, three times the asymmetry of Ref. []T6|| , with opposite sign. This prediction 
is equivalent to a 10 standard deviation effect for the experiment of Ref. [JTJ]. 

While such a large result is anticipated from the estimates given above, it is worthwhile 
to provide an analysis of the amplitudes. These results can be qualitatively understood by 
looking at the interference of the most important low partial wave amplitudes (denoted 
in the spectroscopic notation by l n j or I^lj] the hat indicates CSB amplitudes) 

a (0) _ a ^ -Q) = 4y| J Rep*s 1 cosfl + 2y/2 Re s* u (p 22 - 7|p i + ^21) cos 9. (13) 

This can be compared with the charge-symmetric cross section 

a(9) + a(ir - 9) =2\s u \ 2 + \\p 2 \ 2 +\p 2 \ 2 cos 2 9. (14) 

(The small amplitude 1 S — > p can be neglected in this qualitative comparison.) For the 
case of qm, the largest two CSB amplitudes dominate so that one can clearly see that the 
two interfering terms above add constructively. The other CSB terms, computed earlier, 
have different phases and are therefore subject to destructive interference. 

Note that the total cross section computed from the present model at threshold yields 
an overestimate of the data by about 50% |23|]. We discuss how this could influence our 
present result. The amplitude Su is dominant, so one could arbitrarily reproduce the data 
for the total cross section by dividing the amplitude by y/1.5. This amplitude enters in 
the second (smaller) term of Eq. (|T3"D and in the denominator, so that the resulting value 
of Afb (for all contributions) would increase by a factor of about 1.4. However, the same 
pion rescattering mechanism appears in the both charge symmetric and CSB amplitudes. 
In that case, the error could cancel out in the ratio Aft,. We take an error estimate from 
the geometric mean of unity and 1.4, which gives a theoretical uncertainty of about 20% 

As the isospin- violating amplitudes are small, they are approximately linear in Sm^ — 
5m N /2. Thus we are able to parametrize Aj b as a linear function of Sm N — 6m N /2, at 
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any energy. At E = 279.5 MeV, using the previous results we write 



A lb = - 28% x (l - °™ lSmN _ Wjj (15) 

The predicted value of A/& depends strongly on the value of quantity 5m jy — 5mjv/2. 
Indeed, v4/6 nearly vanishes if this is 1 MeV instead of 4 MeV. However, there are estimates 
cited in the reviews [[I], 0, |TTJ which would obtain 5m n — 5m^j'l significantly greater than 



4 MeV. 

The energy dependence of the amplitudes around this energy is as expected (s wave 
independent of momentum, p wave proportional to q). The asymmetry is not quite as good 
in its energy dependence as the elementary amplitudes: below 300 MeV it is reasonably 
proportional to q; however, above this energy it begins to curve down. With the quark 
model estimate, A fb increases to 1.76% at 300 MeV and 1.87% at 320 MeV. 

The principal conclusion of the present paper, based on Eq. flTSTj is that the successful 
measurement of the value of the forward-backward asymmetry of Eq. (H) for the np — > dir 
reaction will provide a serious test of our understanding of QCD in terms of chiral effective 



field theory. We look forward to the results of Ref. |Tj] with great anticipation. 

We close the paper by suggesting that other experiments probing isospin violation could 
supply additional information, in the whole providing an even more stringent test. We 



have already mentioned the need for more ttN data. Another experiment ||22|| , planned 



at IUCF, would measure the cross section for dd — > air . Charge symmetry prevents the 
reaction dd — > an from occurring, so that any non-zero cross section must be due to the 
breaking of the symmetry. The cross section for dd —>■ an is proportional to the square 
of the matrix element of the CSB pion-production operator. Thus the cross section 
is very small and hard to measure. However, there is no interference with amplitudes 
that respect charge symmetry, which would simplify the interpretation of any measured 
cross section. Furthermore, purely electromagnetic effects are very small, and symmetries 
— parity, angular momentum conservation and the identical nature of the two initial 
deuterons — forbid the production of a pion in a p wave. These features simplify the 
analysis of the reaction dd — > an . 
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